Proton number fluctuation as a signal of the QCD critical end-point 
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We argue that the event-by-event fluctuation of the proton number is a meaningful and promising 
observable for the purpose of detecting the QCD critical end-point in heavy-ion collision experiments. 
The long range fluctuation of the order parameter induces a characteristic correlation between pro- 
tons which can be measured. The proton fluctuation also manifests itself as anomalous enhancement 
of charge fluctuations near the end-point, which might be already seen in existing data. 



The event-by-event fluctuations in heavy-ion collisions 
carry information about the degrees of freedom of the 
created system and their correlations Q]. In particular, 
thermodynamic properties of QCD can be inferred from 
event-by-event fluctuation measurements 0, 0, 0, IE IE ■ 

Of particular interest are fluctuations originating from 
the QCD critical end-point |E II II H EE iH ■ Since 
the fluctuation of the order parameter induces character- 
istic correlations among particles, in particular pions, it 
is expected that the end-point affects the event-by-event 
fluctuations of certain observables in a nontrivial way 
|3,l,i3.J. 

Here we discuss a new observable which may serve as 
a signal of the end-point; the event-by-event fluctuation 
of the net proton number, i.e., the number of the protons 
minus the number of antiprotons observed 0| . 

Our starting point is the fact that the baryon number 
susceptibility xb Qi |l ^ll^ ll^ E, diverges at the 
critical end-point [j, ^ [l^, |20|- Xb is related to the 
average magnitude of the fluctuation SB of the baryon 
number: 



XB 



(1) 



where V and T are the volume and the temperature. 
The divergence of xb is a consequence of the fact that 
the critical point is the end-point of a line of first order 
phase transitions, which are characterized, in particular, 
by a jump in the baryon number density B/V |3lj |. 

If, in a heavy- ion collision experiment, we could mea- 
sure all the baryons, the enhancement of the event-by- 
event fluctuation of the baryon number in a given sub- 
volume would be a signature of the end-point. However, 
about one half of the emitted baryons are undetected neu- 
trons which certainly contribute to the fluctuation of the 
baryon number. To what extent does the proton number 
fluctuation alone reflect the divergence of xs? 

This letter is devoted to clarifying where, in the ob- 
served quantities, the divergence occurs and advocating 
the proton number fluctuation as a sensible and promis- 
ing observable for the search of the critical point in the 
heavy-ion experiments. 



In this work, we confine ourselves to equilibrium ther- 
modynamic fluctuations. Various important issues such 
as the nonequilibrium evolution of the fluctuations will 
be (and some already have been) studied separately. 

For simplicity and clarity we shall work in QCD with 
exact isospin invariance. The relevant corrections due 
to isospin breaking are small as we discuss below. Let 
us first show that in this case the isospin number sus- 
ceptibility, x/i is finite at the end-point. The proof is 
based on the fact that the singular behavior of thermo- 
dynamic quantities near the critical point is due to the 
divergence of a certain correlation length. It is the corre- 
lation length in the cr-channel, the channel with quantum 
numbers of the chiral condensate (V'V') lEIEl- A density- 
density correlator, such as xi — i^/T) J dP'x (Vo(a;)Vb(0)) 
can diverge only if the density can mix with the a field. 
For the isospin density this mixing is strictly forbidden 
by the SU(2)y (isospin) symmetry. The isospin density, 
Vo{x), transforms as a triplet, 3. On the other hand, 
CT is a singlet. The mixing is forbidden and there is no 
singular contribution in xi- [32] 

Small explicit breaking of the SU(2)y symmetry by the 
quark mass difference uiu — md or the isospin chemical po- 
tential /z/ will induce singularity in since m„ — and 
fi] are SU(2)v' triplets [23 j and can produce Vqct mixing. 
In the context of heavy ion collisions corresponding sin- 
gular contributions are negligible. 

We summarize by writing the singular parts of the 
baryon and isospin number susceptibilities: 



XB ~ C 



X/ = (singular parts only). 



(2) 



where ^ is the divergent correlation length of the sigma 
field: ^ — 1/ma-. The zero in eq. (O neglects small 
isospin breaking terms as well as finite terms. The univer- 
sal values of the exponents are given by 7 « 1.2, v « 0.63, 
a w 0.12 < 7. Note that j/iy = 2 - 77, where rj « 0.04. 

Let us discuss the implications of (O. In particular, 
let us consider charge susceptibility. Study of charge 
fluctuations in heavy-ion collisions has attracted much 
attention recently. It was proposed that these fluctua- 
tions might reflect thermodynamic conditions earlier in 
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the collision history, due to charge conservation 0, Q- 
The charge fluctuation magnitude per entropy is a mea- 
sure of the charge per particle or degree of freedom. In 
the quark-gluon plasma (QGP) the charge per degree of 
freedom is smaller. If the charge fluctuation is equilibrat- 
ing too slowly, the observed value will be smaller than 
the equilibrium magnitude which can be calculated us- 
ing the resonan ce g as |^. However, the estimates of the 
charge diffusion [2l| suggest that for the rapidity windows 
achievable in present experiments charge diffusion is very 
efficient in equilibrating charge fluctuations, thus practi- 
cally washing out the "history" effects. Experimentally, 
the QGP suppression of the charge fluctuation is not seen 
[23, 123 , which is consistent with the diffusion estimates 
|2]] |. The effect of the critical fluctuations we are dis- 
cussing here is crucially different from the QGP suppres- 
sion. While the latter is the "history" effect, the critical 
fluctuations are the equilibrium fluctuations pertaining 
to the freezeout point, and the diffusion is necessary to 
establish them. 

The measure of charge fluctuations, the charge number 
susceptibility, %q, can be expressed in terms of xb and 
Xi using the relation Q — B/2 + and the fact that 
isospin symmetry requires {6BSI3) = 0: 



XI- 



(3) 



Eq. Q then implies that the charge susceptibility di- 
verges at the critical point, due to the divergence of xb- 
We now wish to relate the susceptibilities xb, Xi a-nd 
XQ to observable particle number fluctuations. For sim- 
plicity, we shall limit our discussion by considering only 
protons, neutrons and pions. Accounting for other par- 
ticles will not alter our conclusions. In the hadron lan- 
guage, the susceptibilities may be written as 



Xb 
XI 

and 





(5) 



Np for the 



where we introduced notation Np^p = Nj 
net proton number fluctuation, with S denoting event- 
by-event deviation from the equilibrium value. Similar 
notations are used for neutrons and pions. 

Now we concentrate on singular parts of the suscep- 
tibilities and ask a question: what does Q imply for 
the individual particle number fluctuations? It is easy to 
check that the following set of relations between singular 
parts of the particle correlators reproduces the correct 
singular behavior given in |2Jl: 



{SNp^p SNp^p) = (SNn-nSNn-n) = {SNp^p SNn-n) 

{SNp.pSN^+_^~) = {6N„^nSN^+_^-) = 

(singular parts only). (6) 

Some of these equations follow trivially from isospin in- 
variance, but some, for instance, the last equation on the 
first line and that on the second line, require a stronger 
condition. Such relations occur naturally if we attribute 
the divergences to the exchange of a sigma meson, which 
is an isospin singlet. Using eqs. © we obtain 
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— {SNp_p6Np_p) 



= 0, 

1 



yrp {SNp~p^^1p-p/ 



-p6Np 
(singular parts only) 



(7) 



Remarkably, the singular part of the charge fluctuation 
comes from the protons. In other words, had we con- 
sidered only contributions from charged pions in xq, the 
singular parts of tt+tt"*", 7r~7r~, tt+tt" correlators (all are 
singular at the critical point 0) would have canceled each 
other. We see also that the proton number fluctuation 
completely reflects the singularity of the baryon number 
susceptibility, which justifies its use as a sensible probe 
of the QCD critical end-point. 

To provide a simple estimate of how large the net 
proton number fluctuation can become near the critical 
point, we begin by calculating the correlator 



(Sripdnk), 



(8) 



where rip is the net proton number in the momentum bin 
labeled by the value p. In addition to the usual statis- 
tical fluctuation, the correlator ||HJ) receives a contribu- 
tion from the effective interaction with the sigma field 



a, Co 



gaPP, where g is the dimensionless sigma- 



nucleon coupling and P is the Dirac field of a proton. All 
fluctuation observables of the protons can be constructed 
from © 3, 13]. 

Near the critical point, the singular term in ||SJ) is repre- 
sented by a diagram of forward proton-proton scattering. 
A straightforward calculation following gives. 



V{dnpSnk) = 



mlT EpEk 



npil-rip)] 



(singular part only) 



(9) 



where m = 940 MeV is the proton mass, Ep = 
and Up = [exp{{Ep^ fj,B)/T} + while 
is the sigma meson (screening) mass. 

Let us compare the singularity in (jS)) to the singularity 
in The exponent 7/z/ = 2 — ?7in||21lis very close to 2 
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(the anomalous dimension ry of the cr field is small) and 
is equal to 2 in the mean field approximation (77 — 0), 
which is the same as the power of l/m„ in 

In a realistic heavy-ion collision environment finiteness 
of the space-time volume severely prevents m„ from van- 
ishing exactly 0, ^| . The smallest achievable value is 
estimated to be around (3 fm)~^. 

One possible concern is that the rescattering in the fi- 
nal hadronic stage washes out critical point fluctuations. 
In this respect, one should bear in mind that the rescat- 
tering in question includes the exchange of the a quanta, 
which, near the critical point, is the source of the fluctu- 
ations we consider. The critical fluctuations are washed 
out if the final (kinetic) freezeout occurs sufficiently far 
from the critical point. In order to see the effect, one 
should dial control parameters (e.g., reduce the size of 
the ions to raise the freezeout temperature) to bring the 
freezeout closer to the critical point. 

In lO, 5 is taken at zero momentum transfer, i.e., off 
the sigma mass shell. In vacuum, g ~ mZ/jr (/^ — 93 
MeV is the pion decay constant) is quite large ~10. For 
the quantitative estimate below, we assume that g does 
not change appreciably from its vacuum value near the 
chiral phase transition. (See, however, (IJ.) 

First, let us assume that Au-Au collisions at RHIC at 
y/s = 130 GeV froze out in the vicinity of the end-point 
and estimate the effect of Q in terms of the unknown 
mass mo-, which is the measure of the proximity of the 
end-point. At the chemical freezeout, T = 174 MeV and 
/is = 46 MeV Integrating over p and fe, we obtain 
the net proton number fluctuation {{SNp-p)"^) and divide 
it by the sum of proton and antiproton numbers (Np^p) : 



{{SNp^ 



p+pi 



1.0 + 0.062 



RHIC 



2 /200MeV\2 



(10) 



The unity on the r.h.s. is the trivial statistical contri- 
bution. Taking 3 = 10 and m„ = 60 MeV « (3 fm)-i 
|l[l3, we find {{5Np_pf)/{Np+p) « 1.7. 

On the other hand, if the end-point were located at a 
value oi hb of order of a few hundred MeV, as inferred 
from simple model estimates and suggested by the 
recent lattice simulation [13, it is possible that the SPS 
freezout is in the proximity of the critical point. Using 
the freezeout parameters (r,^s)=(168 MeV, 266 MeV) 
at SPS I2I we obtain 



{{SNp-pf 



p+pi 



SPS 



2/200MeV\2 



j ■ (11) 



Note that the coefficient of the second term has a much 
bigger value than in (|10|l . This is because the singular 
term given by grows as the square of the net pro- 
ton number and also because there is a partial cancella- 
tion between protons and antiprotons at RHIC. We stress 
that the main feature in (|10|l and (|ll|l is the singular de- 
pendence on rria, which makes the effect large when the 
freezeout occurs near the critical point. There are other 
effects, which contribute to the r.h.s. of H10|) and (|ll|l . 



but which are not singular near the critical point (e.g., 
initial volume fluctuations caused by impact parameter 
fluctuations). 

Experimentally, separating protons and measuring 
proton fluctuations is a feasible task in the RHIC as well 
as the SPS detectors. We hope that such data analysis 
will be available soon. 

In order to test our ideas on the existing data we 
can, using H10() and Hll(l . estimate the contribution of 
the proton fluctuation to the total charge fluctuation 
characterized by D = A{{5Qf) / {Ntot) H- As Q) 
shows, pions do not contribute to the singular part of 
the charge fluctuation, but they dilute such a contribu- 
tion of the protons. Using {N^^+j^^^-) 10(A^p+p) and 
{{5N^+_^-f)/{N^++^-) « 1-0.3 = 0.7, where the neg- 
ative contribution —0.3 is due to the resonance decays j5|, 
we obtain 



D 

T 



mi 



mp~ 



{Np, 



+ +7r- 



0.: 



(12) 



where we have neglected the cross terms between 5Np and 
5N^± . We see that the fluctuation anomaly in the proton 
sector can result in a larger charge fluctuation than the 
resonance gas value « 0.7 by about 10%. At SPS this 
effect is even stronger. At SPS, using (fTT|l . g — IQ and 
= 200 MeV we get {{5Np_pf)/{Np+p) k. 2.5, and 
with (iVjr++^-) ft! 5{Np+p) we obtain: D/A « 1. 

Before these estimates can be compared to experiment, 
one must take into account the effect of limited accep- 
tance of a given detector. It is easy to see that this effect 
reduces deviations from D — A. Its estimates range from 
few percent corrections ^ to almost complete elimina- 
tion of deviations from D = A [S^l , depending on the as- 
sumptions on the rapidity correlator of fluctuations and 
the width of the acceptance window. We do not discuss 
these issues here, and refer the reader to the literature. 

Experimentally, the data from RHIC suggest that the 
magnitude of the fluctuation is slightly larger than a ther- 
modynamical fluctuation in a resonance gas |22|. This 
effect is even more pronounced at SPS [l^. There are, of 
course, a number of possible explanations, for example, 
(i) acceptance, as we have just discussed [23; (ii) remnant 
initial state correlations; (iii) decay of multiply charge 
clusters; (iv) other nonequilibrium fluctuations (e.g., by 
a mechanism similar to jl^l); In this letter we wish to 
point out that, quite independently of these other ef- 
fects, an equilibrium critical fluctuation due to the prox- 
imity of the end-point could explain the enhancement of 
charge fluctuations observed at RHIC and SPS. The in- 
dependent measurement of proton fluctuations that we 
suggested would be necessary to confirm and sufficient 
to rule out this effect. 

Is it possible that the light sigma effect is seen in both 
RHIC and SPS experiments? If this happens, the region 
of /is where rria is small (< 200 MeV) is rather wide 
- of order hundred MeV. Although unlikely, this might 
not be completely unnatural if one takes into account the 
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fact that VTLa is suppressed on the crossover hne stretch- 
ing from the end-point to fj,B — axis, even though it 
vanishes only at the end-point (see, e.g.. Fig. 5 in 12]). 

As it should be clear from our discussion, measuring 
the charge fluctuations is not the most efhcient way to 
search for the end-point, although the effect may be seen 
in such observables too. A direct measurement of the pro- 
ton number fluctuation as a function of the ^/s of the col- 
lision is both feasible and is less afflicted by other effects. 
Correlation of such a measurement with other proposed 
signatures of the critical point (such as pt fluctuations 
13) would affirm the discovery of the QCD critical point. 

In conclusion, protons carry both the baryon and the 
electric charges. They are sensitive to the fluctuation of 
the order parameter. Due to the peculiar nature of the 
end-point - isospin blindness of the sigma field - the sin- 



gularity of the baryon number susceptibility is completely 
reflected in the proton number fluctuation. Thus the net 
proton number fluctuation is a very useful observable. 
By studying the /i^ dependence of this fluctuation one 
may discover and determine the location of the critical 
point on the phase diagram of QCD. 
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